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Abstract 

Q" 

We introduce a new parameterisation of the effect of unknown corrections 
from new physics on quark and lepton mass matrices. This parameterisation is 
, used in order to study how the hierarchies of quark masses and mixing angles are 

I/-) . modified by random perturbations of the Yukawa matrices. We discuss several 
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examples of flavour relations predicted by different textures, analysing how these 
relations are influenced by the random perturbations. We also comment on 
the unlikely possibility that unknown corrections contribute significantly to the 
hierarchy of masses and mixings. 

43 1 Introduction 

> : 

The observed quark masses and Cabibbo-Kobayashi-Maskawa (CKM) mixing angles 
show a hierarchy m t 3> m c 3> m u , m b ^> m s ^> m d} \V US \ ^> \V c b\ ^> \V u b\ that 
points to an underlying structure in the Yukawa couplings of up and down-type quarks 
to the Higgs boson. This fact has motivated the introduction of several models of 
mass matrices trying to reproduce the experimental values of the quark masses and 
CKM matrix elements using well-structured Yukawa patterns. In general, these models 
consider that at very high energies, and for some reason unknown at present, a weak 
basis is privileged among the infinite set of equivalent weak bases related by rotations 
of the quark fields in flavour space. In this privileged basis, some flavour symmetry 
reduces the number of free parameters in the Yukawa matrices, leading to relations 
among masses and mixing angles that can be experimentally tested. Examples of these 
models include patterns with texture zeros |TJ, democratic Yukawa couplings or 
other symmetries (see for instance 0, 
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At lower scales, where the underlying flavour symmetry is broken, the Yukawa 
textures are modified by renormalisation group (RG) evolution and other radiative 
corrections, with several contributions from Standard Model (SM) and new physics. 
These corrections are unknown, and constitute what is called "anarchy". The "anar- 
chy" may in principle lead to modifications in the predictions of the original patterns. 
The stability of these predictions, when the new contributions are included, can be 
estimated with the addition of random perturbations to the Yukawa matrices in the 
Lagrangian. Of course, the new physics contributions are not random in nature, and 
could be computed if the complete theory was known. At any rate, a simulation with 
random perturbations shows which quantities are stable and which are not, and to what 
extent they are modified in this case. It is already known that the strong hierarchy of 
masses and mixings is not likely to be originated from "anarchy" ||. Recently, it has 
been argued |J that in the quark sector the possible contributions of the latter must 
be very small, in order to preserve a hierarchical structure. However, this result has 
been obtained for a particular texture of Yukawa couplings and in a particular weak 
basis. It would be desirable to obtain a basis and model-independent result indicating, 
at least qualitatively, how the "anarchy" could change the hierarchy of quark masses 
and mixings. 

In the lepton sector the situation is rather different. The masses of the charged 
leptons have been precisely measured, but the knowledge about neutrino masses and 
mixings is still rather poor. After the measurement of atmospheric neutrino oscillations, 
with sin 2 # atm > 0.92, Am^ ta = (1.6 - 3.9) x 1(T 3 GeV [0], it is clear that the pattern 
of lepton masses and mixings is completely different from the one in the quark sector. 
In this direction, it has been suggested that the neutrino mass matrices might be 
generated just from "anarchy" (see also Ref. A complete analysis is not yet 
possible in the lepton sector, for several reasons: (i) some of the mixing angles have not 
still been measured; (ii) the neutrino masses remain unknown and only the differences 
of squared masses can be extracted from oscillation data; (ra) the Dirac or Majorana 
nature of the neutrinos is yet to be determined, and in the first case Majorana-type 
mass terms are not present. 

Our aim here is to study how the "anarchy" would affect a hierarchical pattern, as 
the one found in the quark sector. We first derive the formal structure of the unknown 
corrections to the mass matrices, both for the case of quarks and leptons. We use 
this formalism to investigate how random perturbations modify a given hierarchy of 
quark masses and mixings. Then we discuss the effect in some relations involving 
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masses and CKM matrix elements, which are predicted in several models of Yukawa 
patterns existing in the literature. Finally, we explore the possibility of a significative 
enhancement of the hierarchy of quark masses due to the "anarchy" . 



2 Anarchy in the mass matrices 

The starting point of our argument will be the observation that physical quantities 
cannot depend on the weak basis chosen for the quark fields. In the SM, the mass 
terms of the Lagrangian are written as 

Cm = ~u L v Y u u R - d L v Y d d R + h.c. , (1) 



u 

) 



with v = 174 GeV the vacuum expectation value (VEV) of the Higgs boson, and Y 
Y d the Yukawa matrices, of dimension 3 x 3 in flavour space. In some SM extensions, 
for instance in the minimal supersymmetric Standard Model (MSSM), the masses are 
originated from the VEV's of two Higgs bosons, v u and Vd, with v^+Vj = v 2 . Defining (3 
by tan (3 = v u /va, the mass matrices for up and down quarks are in this case v sin (3 Y u 
and v cos f3 Y d , respectively, with an extra /3-dependent factor. We later comment on 
how this possibility modifies our analysis. 

It is well known that under a change of basis 





, u R = U u R u' R , d R = U d R d' R (2) 

the Yukawa matrices transform as 

Y u -> U{ Y u U% , Y d -> U\ Y d U d R . (3) 
Let us assume that some perturbations are added to the original Yukawa matrices, 

Y u -> Y u + 5Y U , Y d -> Y d + 5Y d . (4) 

The matrices 8Y U , 5Y d are functions of Y u , Y d and other SM and new physics pa- 
rameters. Under the change of basis in Eqs. (|3|), the perturbations must transform 

as 

5Y U ^ul5Y u Ul , 5Y d -> \j{ 5Y d U d R , (5) 
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since the physical observables must be independent of the choice of weak basis Q These 
transformation laws imply that the perturbations have the form 

5Y U = X U Y U + Cu Y u Y u] Y u + r] u Y d Y d ^Y u + ■ ■ • , 

5Y d = X d Y d + ( d Y d Y d] Y d + Vd Y u Y'^Y d + ■■■ . (6) 

The dimensionless coefficients Aj, Q, r]i may be functions of Y u and Y d that are invariant 
under SU(3) flavour rotations, e.g. tr Y U Y U \ tr Y Y d \ etc. and of other SM and new 
physics parameters as well. The higher-order terms are expected to be smaller (if 
perturbation theory is valid), thus in the expansions of Eqs. © the products with five 
or more Yukawa matrices have been omitted. The effect of the \ terms in Eqs. (^j) is 
to rescale the masses by common factors (1 + \ u ) for up quarks and (1 + A^) for down 
quarks, without affecting the hierarchy and the mixing. The Q terms also rescale the 
masses, but with different factors for each quark, m q — > m q (1 + Qm q /v 2 ), and then 
they modify the hierarchies. The rji terms are the lowest-order ones that modify the 
CKM matrix. 

In principle, in SM extensions there may exist a matrix X (not necessarily square) 
of couplings between the quarks and other particles, either transforming on the left or 
on the right side as one of the Yukawa matrices. For instance, if X transforms under 
the change of basis in Eqs. (|3|) as 

X^UlXU$, (7) 

this matrix would originate terms XX^Y U and XX^Y d in Eqs. However, if dia- 
grams involving X couplings gave significant corrections to the Higgs Yukawa vertices 
with a new flavour structure, analogous diagrams with a photon, gluon or Z boson 
would give similar contributions to flavour- changing processes, which are experimen- 
tally very suppressed. Thus we disregard this possibility in the following. 

It is worthwhile noting that the formal structure of the perturbations to the Yukawa 
matrices, derived here from weak-basis independence arguments, coincides with the 
expression of the RG equations in the SM. Setting for instance in Eqs. (^) 

17 o 9 



A " = 16^ U ' 



3Y"Y« + 3Y d Y d ^ + Y e Y^] - (-gj + -g 2 2 + 8gf) } St 



1 These transformation properties for 5Y U and 5Y d do not assume that the Lagrangian is invariant 
under the transformations in Eqs. (^) alone. Within the SM, the Lagrangian is invariant under these 
transformations, but this does not happen in some of its extensions, for instance in the MSSM. Besides, 
at very high energies, some symmetry might single out a special weak basis. Below that scale, and in 
particular at low energies, this symmetry is broken. 
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the SM one-loop RG equations for the Yukawa couplings are recovered |10|, |TT| []. In 
Eqs. (H), Y e is the Yukawa matrix for the charged leptons, g 3 , g 2 and gi are the coupling 
constants of the gauge group SU(3) (g> SU(2) ® U(l) and t = log/i the logarithm of the 
renormalisation scale. 

In the lepton sector the analysis is more involved, because of the possible presence 
of Majorana fermions. In this case the mass terms of the Lagrangian read 



£ m = ~e L v Y e e R - - (u L v R ) 



M L vY 
v Y vT M, 




(9) 



The Dirac mass matrices v Y e and v Y v arise from Yukawa couplings to the Higgs boson. 
The Majorana mass matrix Mr can be included as a bare term in the Lagrangian, 
because the right-handed neutrinos are singlets under the gauge group. The Majorana 
mass term Ml involving the left-handed neutrino fields can only exist at tree-level in 
a renormalisable Lagrangian if a Higgs triplet is present (this is practically excluded 
by precision electroweak data), and Ml is usually set to zero. Both Ml and Mr are 
symmetric matrices. We define Ml = AlMl, Mr = ArMr, with and Ar constants 
with the dimension of mass, in order to express Ml and Mr in terms of dimensionless 
matrices. 

Under a change of weak basis 



(•l 



Vr 



Vr v'r , 



eR 



R °i? 



(10) 



we have 



Y e 
M L 



Vl Y e V* R 



VIMlVI 



Y v - 

Mr 



Vl Y v Vr* , 

> V v r T M r V»r . 



(11) 



2 Note that our Yukawa matrices Y correspond to in Refs. |TT[ . 
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The computation of all the products up to order three transforming as Y e , Y v , Ml and 



M R can be easily done with the programs in Ref. [I2~|. It is found that the perturbations 
to these matrices can be expanded as 

6Y e = X e Y e + Ce Y e Y e ^Y e + r] el Y u Y u ^Y e + 7] e2 M L M L Y e + ■■■ , 
5Y U = X u Y v + ( u Y u Y uJ( Y v + 7] ul Y e Y e ^Y u + rj u2 M L M L Y V 

+U Y V M R M R + U M L Y»*M R + ■ • • , 
5M L = X L M L + C L M L MlM L + r ]Ll (Y e Y e ^M L + M L Y e *Y eT ) 

+VL2 (Y V Y^M L + M L Y v *Y vT ) + xl Y v M R Y vT + • • • , 
5M R = \ R M R + ( R M R MlM R + r ]R (M R Y^Y» + Y» T Y»*M R ) 

+XRY vT MlY v + ••• , (12) 

where the A, (, r], x coefficients are functions of Y e , Y v , Ml, M r and of the 
coupling constants, invariant under the transformations of Eqs. (]TT|) . Under additional 
assumptions, the number of terms can be reduced. For instance, assuming that Ml = 
at tree-level, the terms with ^ e2 , 77^2, £1/2, Xl, (l, Vli, Vl2 and xr can De dropped from 
the expressions. Notice that even setting Ml = in the Lagrangian, our argument 
allows a non-vanishing SMl = XiX v M R Y vT to be generated, though some symmetry 
in the Lagrangian may imply xl — 0. The analysis of the effects of the "anarchy" in 
the lepton sector will be presented elsewhere Il3| . 



3 Effects of anarchy in the hierarchy of masses and 
mixing angles 

We first study how the "anarchy" may change a hierarchical structure in the Yukawa 
couplings. For our discussion we take as benchmark the SM Yukawa matrices in the 
MS scheme at the scale Mz (since the expressions in Eqs. (|]) are basis-independent, 
we can use any basis for the evaluations). The values used for the quark masses are 
collected in Table |I| For the CKM matrix we use, in the standard parameterisation 



I4J, \V US \ = 0.2224, \V ub \ = 0.00362, \V cb \ = 0.0402 and 5 = 1.014. 

We use the perturbations to the Yukawa matrices given in Eqs. @, with random 
coefficients Aj, Q, rji. The moduli of these six independent parameters are generated 
with a Gaussian distribution centred at zero, and for simplicity we assume that the 
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m u 0.0016 

m c 0.68 

m t 175.6 

m d 0.0033 

m s 0.067 
m b 2.9 

Table 1: Quark masses (in GeV) used in the numerical evaluations, 
standard deviations coincide: 

(iAii 2 )' = <ior>* = <w a >'=«. as) 

This is not a serious bias in the analysis, because the moduli of the random parameters 
are not fixed to be all equal, and only the standard deviations of the distributions are 
assumed to be the same. The phases are generated uniformly between and 2tt. The 
only effect of the A, terms is to change the ratio mb/m t ] the mixing angles and the 
ratios of masses of quarks of the same charge are not modified by them. Here it is 
worth remarking that the situation when two different scalars give masses to the up 
and down-type quarks can be reproduced with the substitutions 

yu _^ yu gin p ^ yd _^yd ^ q ^ 



in Eqs. @. This is equivalent to considering 

<|A,| 2 >* 

(IC«I 2 )^ = (M 2 )^ 
(IOI 2 )^ = (\vu\ 2 )^ 



= k sm 2 p , 

= ncos 2 p (15) 



instead of Eqs. ( |13"D . The effect of tan j3 is to modify the standard deviations of the 
distributions of the random parameters. The qualitative behaviour with two scalars 
is the same as in the case discussed below, with only one Higgs. The quantitative 
behaviour is similar, as long as tan/3 ~ 1. 

In the study of the consequences of "anarchy" the following procedure is applied: 
we fix a value of k and generate a set of random matrices, with a number of elements 
between 2000 (for k ~ 0) and 8000 (for k ~ 1). We then select some quantity, for 
instance the ratio m c /m t , and examine its distribution over the set. The la limits 



7 



on this quantity are defined as the boundaries of the 68.7% confidence level central 
interval, evaluated from the sample of random matrices. We begin the discussion with 
this example of m c /m t . We find that this quantity remains fairly stable, even for 
relatively large values of k, e.g. k — 1. In Fig. |] (left) we plot the la limits on m c /m t , 
for k between and 1. In Fig. [I] (right) we plot the distribution of the values of m c /m t 
for k — 1. 




0.25 0.50 0.75 1 0.002 0.004 0.006 0.008 0.010 



K m c /m t (K=l) 

Figure 1: Effect of the random perturbations on the ratio m c /m t . 
Several comments are in order: 

1. The distribution of m c /m t is peaked around the unperturbed value (m c /m t )o = 
3.9 x 10~ 3 , because the random parameters are generated with a Gaussian distri- 
bution centred at zero. It spreads over values larger and smaller than (m c /m t )o, 
but with larger k the smaller values are favoured. This can be observed in both 
plots, and evidences an average increase of m t with respect to m c , due to the 
larger Yukawa coupling of the former. The ratio m c /m t is modified mainly by 
the ( u term, giving corrections 5m c = ( u ml/v 2 , 5m t = ( u m^/v 2 that are much 
larger for the top quark. 

2. The spread of the values of m c /m t increases with k, as it can be expected. The 
growth is linear for small k but it is attenuated for n > 0.5, as can be observed in 
Fig. |I] (left). The upper limit on m c /m t reaches a "saturation" value for k > 0.5 
and does not increase beyond this value. This can be understood as follows: 
As the corrections are larger for the top quark, the ratio m c /m t will generically 
be smaller, and the only chance to have a larger m c /m t is with a cancellation 
between m t and Sm t . For k sufficiently large, the probability of this fine tuning 
to happen is practically constant. 
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3. Most of the values of m c /m t remain close to their original value. Even for k = 1, 
68.7% of the distribution is between the values 2.0 x 10~ 3 and 5.0 x 10 -3 , very 
similar to (m c /mt)o- It is noticeable the long tail of the distribution in Fig. [I] 
(right). For k — 1, the largest value of m c /m t found in the set is 0.15, for which 
the mass hierarchy is washed away. Such values are reached only in an extremely 
small fraction of the sample. 



These plots must not be interpreted as providing any limit on the size of the "anarchy" , 
on the basis of the experimental value of m c /mt- Instead, their meaning is that the 
ratio m c /mt is very stable under perturbations and the mass hierarchy is maintained 
(due to points 2 and 3 above): from an original value (m c /m t ) = 3.9 x 10~ 3 and with 
K = 1 we obtain ratios between 2.0 x 10~ 3 and 5.0 x 10~ 3 most of the time. 

The effect of the random perturbations in the ratio m s /rrib is practically the same, 
as can be observed in Fig. |2|. Although with different numerical values, this ratio shows 
the same behaviour under random perturbations as m c /m t , and the above comments 
apply in this well. The leading correction to m s /mb is due to the rjd term, 

5m s ~ r]dm s m 2 /v 2 , 5mb ~ r]dmbm 2 /v 2 , and the effects are analogous. 




0.02 0.03 0.04 

m./m (k = 1) 



Figure 2: Effect of the random perturbations on the ratio m s /mb. 

The analysis of the two mass ratios m u /m c and rrid/m s shows that they do not 
change when random perturbations are added to the Yukawa matrices, remaining at 
the values (m u /m c )o = 2.3 x 10~ 3 , (md/m s )o = 0.049. This is the same behaviour 
as under RG evolution [jn|, where it is found that these ratios depend weakly on the 
scale. Here we find examples where rrid/m s raises to 0.10, but this only happens in 
a negligibly small fraction of the sample. Therefore, the mass hierarchy between the 
first and second quark generations is also maintained. The remaining independent 
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ratio rrib/mt exhibits a different behaviour, as can be readily noticed in Fig. |3|. The 
maximum of the distribution is displaced to smaller values of m&/m t , and additionally 
there is a long tail for larger values of this quantity. This can be understood as follows: 
Setting \Vtb\ = 1, Eqs. (f|) imply for the third generation 



5mt 
5mb 



A„ m t + ( u rrit/v 2 + r\ u m t m 2 b /v 2 , 
Xdm h + ( d ml/v 2 + r] d m b m 2 /v 2 . 



(16) 



The enhancement for small values of m^/irit is due to the large term £ u mf/t> 2 , while 
the long tail is a consequence of the term rjdmbmf/v 2 , also important. The change of 
this ratio with respect to its original value is limited, though more pronounced than in 
the previous cases. 




0.02 0.04 0.06 

m /m (K = 1) 



Figure 3: Effect of the random perturbations on the ratio m^/mj. 

The modulus of the CKM matrix element V us is not affected by the random pertur- 
bations. As happens for the ratios m u /m c and md/m s , there are fine-tuned values of 
the random parameters for which the corrections yield a much larger value, for instance 
\V US \ ~ 0.4. We stress that this only happens in an insignificant fraction of the sample, 
and for the rest we find \V US \ almost fixed at its original value. The CKM phase S 
does not change either. The moduli of V c b and V u b change when the perturbations are 
added, but their ratio remains constant. (Besides, this is the same behaviour with 
RG evolution that is found for V us , V c b, V u b and the phase 5 |15[| .) The distribution of 
values of \V c b\ can be seen in Fig. |, and coincides with the distribution of \V u b\ up to 
an overall factor. 

Comparing these plots with Figs, [l] and ||] we observe that the effect of the "anarchy" 
on the moduli of the mixing angles V c b and V u b and in the ratios of masses m c /m t , 
m s /mb, m u /mt and md/rrib is similar. This fact is due to the formal structure of the 
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JLfdl 



1 0.02 



0.04 0.06 



Figure 4: Effect of the random perturbations on \V C 



cb\ 



random perturbations, and has non-trivial implications in some flavour relations in the 
next Section. 



4 Effects in flavour relations 



We have selected nine simple flavour relations predicted by several models in the liter- 
ature. Here we do not pretend to test whether these relations are actually fulfilled by 
experimental data. Such analysis ought to be carried out at a very high energy scale, 
of the order of the unification scale Mx ~ 10 16 GeV (where a hypothetical flavour 
symmetry is unbroken), using the Yukawa matrices at that scale. Instead, our aim is 
to investigate whether these relations are modified by the "anarchy". This analysis 
can be safely done at the scale Mz- We prefer to perform the simulations at this scale 
because the extrapolation to the unification scale Mx depends on the particle content 
of the theory between Mz and Mx, and therefore it is different in the SM and its 
extensions. For instance, in the MSSM the running of some parameters, like |V^,|, 
depends strongly on the parameter tan (3. 



Relation I is PL 17 



lKs| = ,p, (17) 

predicted for instance by the five textures proposed in Ref. There are additional 
small corrections to this equality that we ignore. As neither \V US \ nor the ratio md/m s 
change under random perturbations, the accuracy (or inaccuracy) of this relation is 
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not modified. Relation II is 



V. 



ub 



cb 



18) 



It is predicted in textures 1, 2 and 4 of Ref. Q. In this case, the ratios V u b/V c b and 
nT-u/ m c do not change, and this relation also remains unaffected. Relation III is 



m,. 



V ub \ = J — 

V m t 



(19) 



predicted in textures 3 and 5 of Ref. |1[ . In order to investigate this relation, we define 
the ratio 

I / m 

(20) 



Rji 



ub 



m t 



The effects of random perturbations on this quantity can be seen in Fig. |5|. (Recall 
that here i?n is evaluated at the scale Mz- The test of this relation on the Yukawa 
matrices must be done at very high energies, using the RG equations for the SM or the 
SM extension under consideration to evolve the Yukawa couplings.) 




0.5 



1 1,5 
R m (K=l) 



Figure 5: Effect of the random perturbations on the ratio R^, defined in the text. 

The plot in Fig. ||] (left) shows that the perturbations change this ratio from an 
initial value (Rm)o = 0-83 to an interval between 0.55 and 1.56. This implies that the 
corrections to the Yukawa matrices might hide an underlying relation given by Eq. (|i~9|). 
Relation IV is @ 




Correspondingly, we define 
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Figure 6: Effect of the random perturbations on the ratio i?p/, defined in the text. 

This quantity is plotted in Fig. §. Since the ratios m u /m c and \V u b/V c b\ are not modified 
by the random perturbations, the behaviour is similar to R^, up to a global factor. 

Relations V and VI H are very similar and involve only quark masses: 

mu _mc md _m ± 
m c m t ' m s mi, 
The impact of the "anarchy" on them can be investigated with the analysis of the 
ratios 



-Rv = 



m u m t 



Ryi = 



m d m h 



m z c rrig 



(24) 



which are plotted in Figs. [7] and |8], respectively. These quantities are (up to global 
factors) the inverse of the ratios m c /m t and m s /mf, studied above. 




a o.o4 - 



0.5 1.0 

R V (K=1) 



Figure 7: Effect of the random perturbations on the ratio Ry, defined in the text. 



Relation VII |}0|, ^TJ is analogous to Relation I in Eq. 

V td 



v t 



Is 



(25) 
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2 3 4 



Figure 8: Effect of the random perturbations on the ratio Ryi, defined in the text. 



Since the ratio \ V u b/V c b\ and the phase 8 do not change with the perturbations, Relation 
VII remains inaltered as well. Relation VIII is |20 



IK 



cb\ 



V2 



m. 



(26) 



In principle, both sides of this equality change when the random perturbations are 
added to the Yukawa matrices. However, the changes are correlated and the relation 



is not affected. The same happens with Relation IX [|20 



W, 



ub 



l m d m s 
2m 2 b 



(27) 



The non-trivial invariance of these relations under the "anarchy" is due to the formal 
structure adopted for the latter, Eqs. @. 

As remarked at the beginning of this Section, the test of Relations LTX must be 
done at very high energies, where the flavour symmetry originating the Yukawa patterns 
is not broken. In this context, we would like to point out that the test of Relations I, II, 
VII, VIII and IX, as well as other relations invariant under random perturbations, can 
provide a better way to discriminate between different models of Yukawa patterns. On 
the other hand, Relations III, IV, V and VI can be influenced by radiative corrections 
from new physics and hence their observation might be unclear. 



5 Enhancement of the hierarchy 

We want to investigate within our framework to what extent the "anarchy" might 
contribute to (or even generate) the observed hierarchy of quark masses and mixings. 
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The guideline of our analysis is the following: we reduce the hierarchy in the parameters 
of the Lagrangian and then calculate the probability that the hierarchy present in the 
SM is reproduced with the random perturbations. We have studied three examples, 
corresponding to an enhancement of the hierarchy in the up-type quark masses, down- 
type quark masses and CKM mixing angles, respectively. 

1. In the first case, we use the same parameters as in Section |3] but replace the 
top mass by m' t = m t /2. In this situation, the ratio of masses m c /m t is larger 
by a factor of two, (m c /m' t ) = 7.7 x 10~ 3 (that is, the hierarchy is reduced). 
Setting k — 1, we generate a sample of 10000 random Yukawa matrices. A good 
estimate of the likelihood to enhance the hierarchy is given by the probability 
that the ratio m c /m' t after perturbations is smaller or equal than the value m c /m t 
obtained with the true top mass. This probability is only 8.3 x 10 -3 . 

2. In the second example, we use the standard parameters but replace the bottom 
mass by m' b = m^jl. With the same procedure, we find that the probability that 
the ratio m s /m' b after random perturbations is smaller or equal to the SM value 
is 0.12, more than one order of magnitude larger than in the previous case. 

3. In the third example, we use the SM values of the masses, the mixing angle \V US \ 
and the phase 6, multiplying \V u b\ and \V c b\ by two. The probability that the 
ratio |Vd,/V^ s | after random perturbations is equal or smaller than the SM value 
is 0.12. 

With these examples we see that the observed large hierarchy between the second and 
third generations is not an effect of the "anarchy" (as discussed in Ref. ||), though for 
down-type quarks it can receive some enhancements from unknown corrections. The 
two mixing angles V u b, V c b can also be reduced by the corrections. 

6 Summary 

We have introduced a basis- independent parameterisation in order to describe the effect 
of unknown corrections from new physics ("anarchy") in the quark and lepton mass 
matrices. With this parameterisation we have explored the stability of some properties 
of the quark Yukawa matrices against unknown corrections. This has been done with 
the addition of random perturbations to these matrices and a statistical analysis of the 
effect of the perturbations. 
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We have shown that the quark mass hierarchies, namely the ratios m u /m c , m c /m t , 
m d/ m s, fn s / m b and m,b/m t are hardly affected by the random perturbations. Of these 
quantities, m u /m c and md/m s remain constant for most of the values of the random 
parameters. The rest exhibit deviations that in average lead to an enhancement of the 
original hierarchy, but without modifying it significantly. We have also analysed the 
effect in the mixing angles, concluding that neither |V^ S |, nor the ratio |T4&/K*|> nor the 
phase 5 change appreciably with the random perturbations. The mixing angles \V u b\ 
and \ Vcb\ show deviations but still preserving the CKM hierarchy \V US \ 3> \Va\ 3> \V u b\. 
For some fine-tuned values of the random parameters, the strong hierarchy of masses 
and mixing angles is removed, but this only happens in a extremely small subset of 
the sample. We have also found that when the size of the random perturbations (the 
k parameter) is increased, the average effects do not grow linearly but their increase 
rate slows down. 

We have selected nine simple flavour relations among quark masses and CKM mix- 
ing angles, predicted by several models in the literature, discussing the effect of the 
"anarchy" on them. We have identified four relations which are affected by the random 
perturbations. If these relations are fulfilled by the Yukawa matrices, as some theoreti- 
cal models predict [l], 18, [19| , the phenomenological observation may be jeopardised by 
the "anarchy": some flavour properties in the Lagrangian might not be apparent and 
the "anarchy" might blur or hide an underlying flavour relation. On the other hand, 
the remaining five relations discussed are not altered by the "anarchy" and hence their 
analysis could provide a cleaner insight into the structure of the Yukawa matrices. 

Finally, we have used our framework to demonstrate that the possibility that un- 
known corrections give large contributions to the hierarchy of masses and mixings is 
very unlikely, at it has been pointed out before in the literature. We have shown that if, 
for instance, one enhances the ratio m c /m t in the Yukawa matrices in the Lagrangian 
(thus reducing the hierarchy), the probability that the corrections bring it down to the 
SM value is very small. 
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